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Abstract 

We investigate simulations for gauge theories on a Minkowskian 
space-time lattice. We employ stochastic quantization with optimized 
updating using stochastic reweighting or gauge fixing, respectively. 
These procedures do not affect the underlying theory but strongly 
improve the stability properties of the stochastic dynamics, such that 
simulations on larger real-time lattices can be performed. 
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1 Introduction 



First-principles simulations for gauge field theories such as quantum 
chromodynamics (QCD) on a Minkowskian space-time lattice represent one 
of the outstanding aims of current research. Typically, calculations are 
based on a Euclidean formulation, where the time variable is analytically 
continued to imaginary values. By this the quantum theory is mapped 
onto a statistical mechanics problem, which can be simulated by importance 
sampling techniques. In contrast, for real times standard importance 
sampling is not possible because of a non-positive definite probability 
measure. 

Simulations in Minkowskian space-time, however, may be obtained using 
stochastic quantization techniques, which are not based on a probability 
interpretation [H [2]. In Refs. [Sj H] this has been recently used to explore 
the real-time dynamics of an interacting scalar quantum field theory and 
SU{2) gauge field theory in 3 + 1 dimensions. In real-time stochastic 
quantization the quantum ensemble is constructed by a stochastic process in 
an additional "Langevin-time" using the reformulation for the Minkowskian 
path integral |5l |6]: The quantum fields are defined on a physical space- 
time lattice, and the updating employs a Langevin equation with a complex 
driving force in an additional, unphysical "time" direction. Though more 
or less formal proofs of equivalence of the stochastic approach and the path 
integral formulation have been given for Minkowskian space-time, not much 
is known about the general convergence properties and its reliability beyond 
free- field theory or simple models [Gj [7]. Most investigations of complex 
Langevin equations concern simulations in Euclidean space-time with non- 
real actions [HI E] • 

In Ref. |1] real-time stochastic quantization was applied to quantum field 
theory without further optimization. For SU{2) gauge theory no stable 
physical solution of the complex Langevin equation could be observed even 
for small couplings. The physical fixed point was found to be approached 
at intermediate Langevin-times, however, deviations occurred at later times. 
The onset time for deviations could be delayed and physical results extracted, 
if the real-time extent of the lattice was chosen to be sufficiently small on the 
scale of the inverse temperature. This procedure provided severe restrictions 
for actual applications of the method. In contrast, for self-interacting scalar 
field theory stable physical solutions were observed. 

In this paper we investigate real-time stochastic quantization for gauge 
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theories employing an optimized updating procedure for the Langevin 
process. We consider optimized updating using stochastic reweighting or 
gauge fixing, respectively. These procedures do not affect the underlying 
theory but strongly improve the stability properties of the stochastic 
dynamics. For SU{2) gauge theory in 3 + 1 dimensions we demonstrate 
that gauge fixing leads already to stable physical solution for not too small 
(3 ~ l/f?^, where large (3 correspond to going to the continuum limit of the 
lattice gauge theory. Where applicable, the results are shown to accurately 
reproduce alternative calculations in Euclidean space-time. In order to gain 
analytical understanding and to compare to exact results we also investigate 
U{1) and SU{2) one-plaquette models. 

The paper is organized as follows. In Sec. [2] we briefiy review real- 
time stochastic quantization for non-Abelian lattice gauge theory following 
Ref. [4j. The U{1) one-plaquette model of Sec. 13.11 is used to introduce the 
concept of stochastic reweighting in Sec. 13.1.31 The simplicity of the model 
allows us to compare simulation with analytical results and to investigate in 
some detail the fixed point structure and convergence properties is Sees. 13.1.41 
and 13.1.51 In Sec. 13.21 we consider the SU{2) one-plaquette model and 
introduce some important notions that will be employed for the optimized 
updating using gauge fixing for the lattice field theory in Sec. HI We present 
conclusions in Sec. [5] and an appendix provides some mathematical details. 



2 Real-time gauge theory 

Gauge theories on a lattice are formulated in terms of the parallel transporter 
Ux,^ associated with the link from the neighboring lattice point x + fiio the 
point X in the direction of the lattice axis /i = 0,1,2,3. The link variable 
Ux,^l = U^^^ _^ is an element of the gauge group G. For G = SU (N) or U{1) 
one has U~j^ = ^, however, since we will consider a more general group 
space in the context of stochastic quantization this will not be assumed. 
Therefore, we keep U~j^^ in the definition of the action, which is described in 
terms of the gauge invariant plaquette variable 

where 11")^^ = Ux,^u ■ The action on a real-time lattice reads 

S[U] = -/3oEE{^(Trt^^,o. + Tr[/-o\)-l} 
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E E { ^ f^-^^- + - 1 } 



(2) 



X «J 



with spatial indices i,j = 1, 2, 3. Here the relative sign between the time-like 
and the space-like plaquette terms reflects the Minkowskian metric, and 

with the anisotropy parameter 7 = ag/at on a lattice of size (Ngas)^ x Ntat- 
Because of the anisotropic lattice we have introduced the anisotropic bare 
couplings go for the time-like plaquettes and gs for the space-like plaquettes. 

Using stochastic quantization the real-time quantum configurations in 
3+1 dimensions are constructed by a stochastic process in an additional (5th) 
Langevin-time O El |2]. For a discretization with stepsize e the Langevin- 
time after n steps is -^n = The discretized Langevin equation for the link 
variable reads with the notation U'^^^ = Ux^^i^&n+i) and U^^^ = U^^^i^n) [1] 



(4) 



where differentiation in group space is defined by 



D^^afiU^,,) = -^f (e*"""t/.,^) U=o (5) 

with the generators of the Lie algebra and a = 1, . . . , N'^ — 1 for SU{N). 
For the action ([2]) one has 
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where we have defined Pij = Ps, Poi = Ao = —Po and 

— TT - TT~^ TI^'^ -I- TI~^ TT^^ TT ~ 
Cx,^u = Ux,vUx+i>,fJJ^j^^^^ + U^_0 j^Ux-u,fiUx+fi~C',i' • (7) 



With Ux,^iCx,^iu = Ux,f,u + Ux,fii-u) and Cx,^uUx}, = U^ ^^ + f/^j^ 



x,fi{—u) 



one 



observes that the sum in Eq. is over all possible plaquettes containing 
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Ux,^ ■ Following Ref. [Ij the Gaussian noise rj^^j^a = Tjx^ai'&n) appearing in (jlj) 
is taken to be real and satisfied 

{Vxi^a) = , {T]x^a rjyub) = 2 ^^^.^xj/^ab • (8) 

Expectation values for observables can be obtained from solving equation 
(jlj) for sufficiently large Langevin-time by performing noise averages or, 
alternatively, from Langevin-time averages |10j . 

For instance, specifying to SU{2) gauge theory \a = aa {a = 1, 2, 3) 
represent the Pauli matrices, and one can make further simplifications 
using Tr(t/-V") = -Tr(?7a") for any element U G SU{2). The latter 
simplification also holds for U G SL{2,C). This is relevant since possible 
solutions of Eq. (jl]) may respect this enlarged symmetry group. Taking 

U^^^ = e'^^^^^"^/^ (9) 

the vector fields A^^a need not to be real for U G SL{2,C). The complex 
matrix A'^^aa still remains traceless, however, the Hermiticity properties are 
lost. As a consequence, it is no longer possible to identify W with as 
is taken into account in Eq. ([2]). This corresponds to an extension of the 
original SU{2) manifold to SL{2, C) for the Langevin dynamics. Only after 
taking noise or Langevin-time averages, respectively, the expectation values 
of the original SU{2) gauge theory are to be recovered. Accordingly, if the 
Langevin fiow converges to a fixed point solution of Eq. (jl]) it automatically 
fulfills the infinite hierarchy of Dyson-Schwinger identities of the original 
theory [Hill]. 

3 Optimized updating: simple examples 

3.1 One-plaquette model with U{1) symmetry 
3.1.1 Direct integration 

As a first example we consider the one-plaquette model with U{1) symmetry. 
For U = e*"^ the action is given by 

^0 = ^{U + U-^) = Pcosif (10) 

^It was suggested in earlier literature to replace Sf^i, on the right-hand side of Eq. ([8]) 
by g^i,. However, in this case solutions of the Langevin evolution would not respect the 
Dyson-Schwinger identities of the underlying quantum field theory, as is shown in Ref. [4] . 
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with real coupling parameter (3. The one-plaquette "partition function" is 

Zo = / d</.e*^« = 27rJo(/5), (11) 

where Jo(/3) denotes a Bessel function of the first kind [12j, with 
Jo(l) ~ 0.765. The average of an observable 0{ip) is obtained as 

(0)o = ^ r^d^e'^°0(^). (12) 

For real (3 the integrand in Eq. (JT2l) is not positive definite, which mimics 
certain aspects of more complicated theories in Minkowskian space-time that 
will be considered below. In contrast to those more realistic theories, the 
one-plaquette model has the advantage that the elementary integrals can 
be performed and the results directly compared to those obtained from 
stochastic methods. 



3.1.2 Complex Langevin equation 

In principle, adding a Langevin-time dependence (p ^p{^n) all observables 
can be computed from a solution of the discretized Langevin equation 

ip = + V^/er] 

= ip — ie P simp + ^/e r] , (13) 
using a notation as in Eq. (jlj). The real Gaussian noise fulfills 

(r^) = , {r^v) = 2 (14) 

according to Eq. ([8]). 

In view of the aim to compute plaquette averages, we consider the average 
of the function e'^'"'^ with integer /|§ We first compute averages by analytic 
or direct numerical integration according to Eq. (fT2l) . and then compare to 
the result for the same quantity obtained from a stochastic process using Eq. 
f fT3|) . For (3 = 1 and / = 1 one obtains from Eq. f|T2|) 

{e^o = ^ ^0.575, (15) 



^We will not investigate here the question of defining roots of group elements from 
stochastic processes and do not consider non-integer I. 
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which can be compared to the result from the solution of the Langevin 
equation (IT^ : 



(e*'^)o °" -0.009(±0.006) -i0.00006(±0.00007). (16) 

This result was obtained for /3 = 1 using a Langevin stepsize e = 10^^ 
from a Langevin-time average over 10^° steps. The error in brackets gives 
the statistical fluctuation of the average. One observes that the simulation 
yields a wrong result that is compatible with zero, in contrast to the non- 
vanishing imaginary value (1151) obtained analytically. A similar failure of the 
stochastic method can be observed for averages of other functions as well. 

3.1.3 Optimized updating by stochastic reweighting 

The very same average values as above may be computed with the help of 
the partition function for a different action Sa , 

27r 

dipe'^-, (17) 







from 

{0)a = ^ rdve^'-0{^) (18) 

Jo 

using standard reweighting techniques. If we define 

then the expectation value f[T^ can be identically written as 

ff"dvpe-^-^„0(<^) {uj^O)^ 

~ ^ ~ in \ ■ ^ ^ 

We will consider the family of actions 

Sa = Sq + a(p = P cos(y9 + aip (21) 
with integeJ^a, such that the reweighting function (IT^ reads 

= e"*"^ . (22) 
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—6.66 


-6.66(±0.03) 





0.0338(±0.06) 



Table 1: Averages (e*''^)^ for different values of the reweighting parameter a and 
fixed 13=1. Listed are results for the real and imaginary part of (6*'*^)^ denoted as 
Re and Im, respectively. Compared are exact results from direct integration with 
simulation results using a stochastic process. For the exact values three significant 
digits are given. For the stochastic method the given error in brackets reflects 
statistical fluctuations. 



For the action ( ]2T1) the discretized Langevin equation is given by 

ip' = (f — ie P simp + iea + y/er] . (23) 

As a consequence, for a 7^ the average value in Eq. (!20l) is computed with 
the help of a different stochastic process than in Sec. 13.11 For (3 = 1 and the 
very same Langevin parameters as above we obtain for a = 1 the simulation 
result: 

(e^'^)o = ^^^°=^ = -0.001(±0.003) +i0.575(±0.0004), (24) 

which is close to the exact result given in Eq. (fTSll . in contrast to the failure 
of the method without reweighting. 

In Table [1] we list results for averages using various values of a and fixed 
(3 = 1. Shown are the results for Re(e*''^)Q, and Im(e*''^)Q, both from direct 
integration ("exact") of Eq. f|T8l) and using a stochastic process ("stochastic") 
according to Eq. fl23l) . where we use the same Langevin parameters for the 
numerics as before. The first two rows correspond to a reweighting parameter 
a = 0, i.e. no reweighting, showing the strong disagreement of simulation and 
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P 

Figure 1: The real and imaginary part of the average {e^^)a=i as a function of (3. 
The hnes represent averages obtained from direct integration, while the symbols 
are measurements using a stochastic process. 

exact results in this case. We obtain similarly bad results as long as a < /3. In 
contrast, Table [T] shows for a = 1 accurate values obtained from simulation. 
We note that in this case the value for (3 agrees with the value chosen for 
the reweighting parameter a, i.e. a = (3 = 1. For a = 2 the results are still 
accurate, which we observe also for several other a and (3 as long as a > /3. 
We typically need to collect substantially more statistics for negative I to 
keep statistical errors small as compared to the case with positive /, which 
will be addressed in Sec. 13.1.51 

The dependence of the accuracy of the simulation outcome on the relative 
size of a and (3 is further illustrated in Fig. [1], where we consider results for 
various values of (3 and fixed a = 1 with / = 1. Shown are the values for 
the averages Re(e*'^)a=i and Im(e*'^)a=i again from direct integration (lines), 
which yields 

as well as using a stochastic process (symbols). One observes that for /5 < a 
the simulation method gives accurate results for this quantity, while for 
somewhat larger (3 they can deviate substantially. In Fig. [2] we show results 
for (e**^) as a function of integer values of a for the case (3 = a, which shows 
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Figure 2: Shown are the real and the imaginary part of the average (e**^) as 
a function of integer a with f3 = a. 



very good agreement between stochastic averages and direct integrations.'^ 
The above analysis shows that the accuracy of the simulation method 
depends strongly on the employed stochastic process. It suggests that for 
a given model there is an optimized stochastic process, which may yield 
accurate results. The conditions that have to be met in order to obtain 
quantitative estimates will be further explained in the following. 

3.1.4 Fixed point structure 

The stationary solutions of the noise averaged Langevin equation (12^ are 
determined by 

- /3{sm(p)a + a = 0. (26) 

Taking into account that the dynamical variable can become complex with 
= fR + if I this is equivalent to the statement that the equations 

{cosipR smhipi)a = 0, 

(sin cosh = -, (27) 

have to be simultaneously satisfied. Neglecting fluctuations, i.e. disregarding 
for a moment the noise averages, Eqs. (l27l) correspond to the classical fixed 
point condition dS/d(p {(p = (p*) = 0. The first of these equations constrains 
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either = yj} = or Lp^ = (f*j^ = tt /2 or 37r/2 for the real part of 
the fixed point value. Taking into account the second equation of ([27]), 
there are two distinct real solutions for j3 > \a\. There are two complex 
conjugated solutions for /3 < |a| and for /3 = |a| there is a single real 
solution. We emphasize that these last statements about the possible fixed 
points just consider the first derivative of the classical action with respect to 
the dynamical variable, instead of the noise averaged quantity {dS/dip). 

The classical approximation becomes exact for /3 = |q;| — > cx). In this limit 
the integral (!20|) is dominated by the value where the oscillatory integrand 
shows slowest variation in ip, i.e. for ip = ip*. As a consequence one obtains 

,^_^£%l^e~^^ (28) 

with (fi* = 7t/2 for the positive sign in the exponent of the integrand, 
i.e. a = P > 0, and ip* = 37r/2 for the negative sign corresponding to 
a = —[3 < 0. We will see below that important aspects of the results for 
finite (3 of Sec. 13.1.31 can be understood already from the classical fixed point 
condition, i.e. neglecting fluctuations. 

For Fig. [3], we have evaluated dS/ dip for the range of values < y^ij < 27r 
and — vr < v?/ < tt and plotted its real and imaginary part as a vector with 
origin at each y)- value. The size of the real and imaginary part of dS/dip 
determines the direction and angle of each vector, however, we normalized 
their length for better visibility. The left figure employs a = 0, for which one 
infers from Eqs. (1271) two classical fixed points at (p* = and (p* = it. From 
the arrows it can be seen that these do not correspond to attractive fixed 
points, where the drift term in the Langevin equation f|T3l) would tend to 
focus the Langevin flow. Instead they are "circular" with opposite rotation 
directions for the two points. Taking into account the 27r-periodicity of 
the dynamical variable one observes that the fixed points are equidistantly 
separated along the real axis in this case. 

These properties of the classical fixed point determine to a large extend 
the full Langevin flow, i.e. the behavior of the dynamical variable in the 
presence of fluctuations due to the noise term ~ in Eq. (123|) . In order to 
visualize the distribution of ip, we make snapshots of the Langevin process 
with a time-step of = 1, and plot the values in the complex plane. The 
resulting distribution for 1 = f3 > a = 0, i.e. without reweighting, is given 
in the left graph of Fig. [31 One observes that the values are rather evenly 
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Re((p) Re((p) 

Figure 3: Shown is the real and imaginary part of dS/dip plotted as a vector 
with origin at each ip-value and with normalized length for better visibility. Here 
we consider /3 > a, for which the stochastic method fails. For the left graph 
1 = P > a = 0, while the right graph employs 1.5 = f3 > a = 1. From the arrows 
one can infer the two non-attractive, classical fixed points on the real axis. Also 
shown is the distribution of f as obtained from the full solution of the respective 
complex Langevin equation. One observes relatively wide distributions with values 
having positive as well as negative imaginary parts. This has to be compared to 
Fig. below for the case P ^ a. 

distributed along the real axis and very loosely centered around it in the 
complex plane. 

For comparison the right graph of Fig. [3] shows the corresponding results 
for 1.5 = j3 > a = 1. One observes the two circular fixed points, which are, 
however, no longer equidistantly distributed along the real axis. Accordingly, 
the distribution obtained from the full Langevin dynamics varies considerably 
along the real axis, with a larger density of points where the classical fixed 
points are closest to each other. The larger value for (3 compared to the one 
employed for the left graph leads to a smaller width of the distribution in 
the complex plane. However, in contrast to the case a = 0, one observes 
that the values for the reweighted theory are predominantly localized in the 
positive-*/?/ half-plane. 

The latter tendency continues when j3 is decreased with respect to 
The left graph of Fig. H] shows results for (3 = a = 1, for which only one 
classical fixed point at a real value appears. This stationary point is neither 
attractive nor repulsive, with all arrows pointing towards the fixed point 
for positive imaginary part of (f and away from it for negative imaginary 
part. From the distribution in the left graph of Fig. H] one observes that 
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Figure 4: Same as in Fig. [3] but for /3 = a = 1 (left) and 0.5 = /? < a = 1 (right). 
For P = a one observes the real fixed point with all arrows pointing towards it 
for positive imaginary part and away from it for negative imaginary part of the 
dynamical variable. For /3 < a an attractive fixed point in the positive-^?/ half- 
plane and a repulsive one in the lower half-plane appears. The comparably narrow 
distribution for ip shows that the Langevin flow spends practically all the time 
near the attractive side with positive imaginary parts. For the considered cases, 
where /3 < a, the stochastic method is found to give quantitative results. 

the dynamical variable spends most of the Langevin-time near the attractive 
side of the fixed point. Finally, for 0.5 = (3 < a = 1 the two complex fixed 
points are visible from the right graph of Fig. HI The one in the positive- 
(pj half-plane is attractive with all arrows pointing towards the fixed point, 
while the other in the negative-yj/ half-plane is repulsive. Accordingly, the 
distribution shown in the right graph of Fig. H] is practically entirely localized 
in the positive-yj/ half-plane. 

Fig. [5] shows the same as the left graph of Fig. H] but with larger 
j3 = a = 10 (left) and j3 = a = 100 (right). As the value for P is increased the 
distribution is more and more centered near the vicinity of the classical fixed 
point. This reflects the fact that the (a = /3)-reweighted one-plaquette model 
has a well defined limit /? — > oo described by Eq. (!28l) . in which fluctuations 
are suppressed. This qualitative property of a suppression of fluctuations for 
large (3 will also be encountered in the discussion for the non-Abelian field 
theory in Sec. HI 
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Figure 5: Same as for the left graph of Fig. H] but for P = a = 10 (left) and 
P = a = WO (right). 

3.1.5 Convergence 

In view of the above findings we consider in the following analytical 
arguments under which conditions convergence of the stochastic process to 
accurate results may be expected. For this we associate the stochastic process 
f l23|) to a Langevin-time dependent distribution Pa{ip]'dn) for the stochastic 
variable ip. Using a notation as in Eq. fl23|) by writing P^iv) = Pa{'^] ^n+i) 
and Pa{f) = Pai,^] ^n) its evolution can be obtained from 



where the brackets indicate noise average according to Eq. (IT^ . Expanding 
the ^-functions and keeping only terms up to order e gives the Fokker-Planck 
equation 




— if — le 




(29) 




(30) 



In the continuum limit e 



we write 



d'd 



HMv)Pai¥^;^), 



(31) 



with the Fokker-Planck " Hamiltonian' 



+ i 




(32) 
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We first consider the limit of large j3 = \a\, where we have seen in 
Sec. 13.1.41 that the stochastic process converges and is properly governed 
by the classical fixed point dS/dip {(p = = 0. According to Eq. fl28l) the 
stationary distribution 



\imP^{v;^) = P:{^) (33) 

17— >00 



in this case is described by 



P:{v) ^ \un — %-^*). (34) 

p=|a|^oo Z/Q, 

In order to study the Langevin-time dependence of averages of an observable 
0{ip) in the limit of large jS = \a\, we consider the difference with respect to 
the stationary solution, i.e. 

A0{^;^) = J dipO{^)AP^{^;^), (35) 

where APq, = P^ — P* for properly normalized distributions. Since for the 
classical fixed point the limiting distribution (I51|) has a compact support 
away from the boundaries of integration, and assuming analyticitj{§ in Lp, we 
may use partial integration to write: 

- J d^o{ip)HM^)APa{v;^) = 

Since we are interested in plaquette averages, we consider again 0{ip) = e*'''' 
with integer I and Sa{<f) given by Eq. (!2T|) . According to the Fokker-Planck 
equation (13T!) the Langevin-time evolution for the observable average (135|) is 
then described by 

^ AO{ip; i9) = - (/2 + la) AO{ip; ^) + l(3 J dip sm{ip)0{ip)AP^{ip; ^) 

~ - {l^ + l[a- f3sm{-^)]) AO{ip;^) . (37) 

For the last approximate relation we used that with Eq. fl34l) the integrand 
is peaked near the classical fixed point such that an appropriate constant Ip 



^See, however, the discussion in Ref. 
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may be found with ^ ~ in order to simplify the remaining integral and to 
get a closed equation for A0(v9; {}). Eq. (1571) has to be evaluated for /? = \a\ 
but in the notation we keep f3 and a separately for further discussion. The 
Langevin-time dependence of the observable average is then given by 

With 

sin(^) = sin(^j:j) cosh(^j) + i cos(^j:j) sinh(^j) (39) 

one obtains from Eq. (138!) a convergent result if 

P + l[a- (3 sin(^^) cosh(^^)] > . (40) 

For P = \a\ ^ oo we have Tp = ip* a.t sufficiently large ^. At the classical 
fixed point dS/dip = a — sm{ip*) = which leads to lim^^oo AO{ip; {}) ~ 
lim^^oo e~' = 0. Accordingly, the stochastic method is expected to 
converge well to the stationary solution in this case, which we indeed 
observe from the full solution of the Langevin equation as described above 
in Sections 13.1.31 and 13.1.41 

We note that the solution fl55]) happens to reflect important qualitative 
properties of the above discussed results also for finite P and for /9 7^ The 
partial integration leading to Eq. (!36|) is based on analyticity arguments and 
a compact distribution for the stochastic variable away from the boundaries 
of integration. The latter is also assumed for the step to the second line of 
Eq. ( 1371) . Of course, in case the dynamics is governed by a classical fixed 
point, i.e. if the distribution of the dynamical variable is centered around 
Tp with dS/dip{ip = ■^) ~ 0, the condition ( HOl) is automatically fulfilled for 
any /3 or a because of the fixed point condition (|27|) . However, fluctuations 
often play an important role and to analytically argue that the Langevin flow 
not only converges but converges to the correct value is more involved if the 
dynamics is not governed by a classical fixed point. 

The importance of fluctuations can be observed, e.g., from Fig. [21 where 
foi' 1^ ^ (9(10) substantial deviations from the classical fixed point value 
occur. If dS/dip{ip = ^) 7^ then the condition (l40l) can be formally written 
as 

I + a 

sin(^^) cosh(^^) < — — (41) 
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for 1^0. This suggests that convergent results might be difficult to obtain 
for (3 ^ a, which is indeed in accordance with our findings of Sec. I3.1.3[ For 
instance, in Fig. [1] one observes accurate results for I = 1 from the stochastic 
method if /3 < a = 1, but a failure of the method for larger j3. This coincides 
with the fact that for /5 < a we find a rather compact distribution for 
as exemplified in Fig. |H In this case one observes that the condition fHTl) 
is approximately verified if Tp^ and Tpj are allowed to take on all (/^-values 
of significant support. In contrast, for (3 larger than a we find rather wide 
distributions as exemplified in Fig. [3] and the assumptions leading to (HT!) may 
not be justified. In these cases it turns out that the observed distributions are 
also difficult to reconcile with condition fl4ip . In particular, the case a = 
with (3=1 (left graph of Fig. [3]) for / = — 1, which was employed in Sec. 13.1.21 
to demonstrate the failure of the stochastic method, is clearly violating that 
condition. We also verified some further details, e.g., that for a = 1 and 
a = 2 the observable averages with / = ±1 and / = ±2 converge faster up to 
a prescribed statistical error if / is positive, as suggested by condition ( |4T1) . 



3.2 SU{2) one-plaquette model 
3.2.1 Direct integration 

We consider a theory where the action with real coupling parameter f3, 

^(f/) = ^Trf/, (42) 

is invariant under the symmetry transformation 

U W-^UW (43) 

with U,W ^ SU (2). For the analytical calculations we use a parametrization 
of SU (2) matrices in terms of an angle and a unit vector n: 

U ((^, n) = e*'^"-"/2 = (^cos |) 1 + ^ (sin ^) n ■ a , (44) 

where < (/? < 27r with the three Pauli matrices a. Using the Haar measure 
in terms of these variables we obtain averages of an observable 0{U) by direct 
integration from 



2 



I / (sin|)"e^/^-^(-/^)0(f/(^,n)) , (45) 
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where Q{n) is the uniform measure on the unit sphere and 



Z = (sin|)%'/^™^(-/2) = jJm- (46) 

For instance, with Tr f//2 = cos((y9/2) the plaquette average as a function of 
13 is 

¥") - 'Hi- 

Note that for this particular observable the value of the integral coincides 
with the one of the plaquette average in Eq. ( 125|) for the (a = l)-reweighted 
U{1) model, and the analytic results are plotted as solid curves in Fig. [H 



3.2.2 Complex Langevin equation 

We will compare results obtained from direct integration to estimates from 
a stochastic process. As described in Sec. [21 this requires an extension of the 
original SU{2) manifold to SL{2, C) for the Langevin dynamics. Writing 

TT 1 , -r - f a + ibs b2 + ibi \ , . 

U = al + ib-a={ , , (48) 
\^ -62 + ibi a - ib^ J 

with det U = a? + b'^ = 1 the coefficients a and b are complex numbers for 
U G SL{2, C). The Langevin equation follows from Eq. (j4]). We expand the 
exponential in that equation to first order in e, which means we must include 
the square of the noise term, which is proportional to unity. The evolution 
equation then reads for the one-plaquette model with action fl42p : 

U' = (^dl + zJ2a,(^-ef^Ti{aM) + V~^VayjU, (49) 

with Gaussian noise rja corresponding to Eq. ([8]). In order to 
stay in group space the constant d in this equation is calculated 
from d = a/1 — (— e/9Tr(craf/)/2 + ^/erjay. Alternatively, one can also 
appropriately normalize that matrix. In contrast to the procedure of Sec. 13. 
Eq. (149|) describes the stochastic dynamics directly in terms of group elements 
U. This is closer to what we will do for the SU{2) gauge theory in 3 + 1 
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Figure 6: The real and imaginary part of the plaquette variable TrU from 
snapshots with constant Langevin-time stepping for the SU (2) one-plaquette 
model with (3 = 1. The left graph shows the wide distribution of values obtained 
from the standard Langevin dynamics, while the right graph displays the compact 
distribution from the (" gauge- fixed" ) optimized process (see text for explanation) . 

dimensions in Sec. IH and is used to introduce some concepts that will be 
employed for the field theory as well. 

Applying the stochastic process to a computation of the plaquette average 
for (3 = 1 yields a result consistent with zero: 



which disagrees with the non-vanishing exact value iJ2(l)/Ji(l) — i 0.261 
given by Eq. fHTl) . The situation is analogous to the one described in Sec. 13.11 
In particular, the distribution for the stochastic variable obtained from the 
solution of the Langevin equation fH9l) exhibits similar qualitative patterns. 
For instance, the left graph of Fig. [H] shows for /3 = 1 the distribution of Tr U 
in the complex plane. One observes a wide distribution, reminiscent of the 
left graph shown in Fig. [3] for the U{1) model without optimized updating. 

3.2.3 Optimized updating by "gauge fixing" 

An optimized updating scheme for the Langevin dynamics with the aim to 
calculate accurate results may be achieved in various ways. Motivated by 
the results of Sec. 13. H the optimized updating should control the growth 
of fluctuations for the complex Langevin equation. Since the SU{2) one- 
plaquette model has a global symmetry described in Eq. fH51) . which is 




without 
optimization 




(50) 
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observable 


exact Re 


stochastic Re 


exact Im 


stochastic Im 


(Trf//2) 





-0.004(±0.006) 


0.261 


0.260(±0.001) 


((W/2)2) 


0.216 


0.217(±0.003) 





-0.001(±0.002) 



Table 2: Results from the optimized stochastic process, as described in Sec. 13.2.31 
are compared to exact values from direct integration for (5 = 1. For the latter only 
three significant digits are given. Listed are results for the real and imaginary part 
of (TrU/2) and ((TrU/2)^) denoted as Re and Im, respectively. 

reminiscent of a local gauge transformation in the corresponding field theory, 
one may use this symmetry to "gauge- fix" certain variables in order to 
constrain the growth of fluctuations. For the complex Langevin equation 
the plaquette variable U as well as the matrix W in Eq. are elements 
of SL{2,C). In the following we will use it in order to diagonalize U 
after each successive Langevin-time step. For the representation (jUj) 
this corresponds to the "gauge-condition" n = (0,0,1), or, with the 
parameterization of Eq. ( l48l) we write after each Langevin updating step 
U' = diag (a -|- za/1 — a^, a — — cl^) , which corresponds to choosing 

b= (o,o,v/r^)0 

Using this procedure we find that, in contrast to the wrong results 
displayed in Eq. (!50|) for /? = 1, the optimized stochastic process now 
reproduces correct averages. This is exemplified in Table [2] for two different 
observables, {TtU/2) and ((Tr[//2)^). One can compare the distribution 
of values obtained from snapshots at equidistant Langevin-time steps with 
and without optimized updating. The right graph of Fig. [6] shows that the 
"gauge- fixing" leads to a compact distribution, in contrast to the one from 
the unmodified process displayed on the left of that figure. Moreover, for 
(Tr[//2) as a function of (3 we obtain a very similar plot than the one shown 
in Fig. [T] with accurate results for /5 < 1 and wrong results for somewhat 
larger /?. Of course, according to Eqs. (!25l) and (H7I) . the expectation values 
have to coincide in this particular case in the range where correct values are 
obtained. 

We note that the discussion can be performed along very similar lines 
than what has been done in Sec. 13.11 for the U{1) one-plaquette model. This 

^Alternatively, one can choose the negative sign with b ~ (0, 0, — Vl — a^), which makes 
no difference for the following discussion. 
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can be observed from the fact that the above optimized updating procedure 
can be mapped to a Langevin process with one degree of freedom only, given 
by 

ip = - «e — sm — + e cot — + V e (,51) 

with white noise 77. This Langevin equation follows directly from Eq. (l46ll by 
writing the reduced Haar measure as an exponential. The equivalence with 
the optimized updating is verified explicitly in an expansion of Eq. (H9l) to 
order e in the appendix. 

In contrast to the (a = /3)-reweighted one-plaquette model of Sec. 13. ![ 
the SU{2) one-plaquette model has no well-defined limit (3 —>■ 00, in which 
fluctuations are suppressed. This can be observed, for instance, from the 
integral in Eq. ( H7|) . For the following discussion it will be an important 
property that the SU{2) gauge theory has a well-defined limit f3 — >■ 00, which 
corresponds to the continuum limit of the lattice theory. 



4 SU (2) gauge theory in 3 + 1 dimensions 

4.1 Stochastic dynamics without optimized updating 

Real-time simulations for quantum field theories at non-zero temperature 
require a complex time-path, where the imaginary-time extent is given by 
the inverse temperature ~ 1/T and the physical time determines the real- 
time extent. For scalar and SU{2) gauge theory this has been investigated in 
Ref. P] to which we refer for further details. Here we re-consider SU (2) gauge 
theory. In this case the real-time quantum dynamics in 3+1 dimensions is 
obtained from a stochastic process in the additional (5th) Langevin-time 
according to Eq. (jlj) explained in Sec. [2l The numerical results will be 
obtained using equal couplings for the time-like and space-like plaquettes 
in Eq. ([31), i.e. gl = gl = g'^ r-. 1/ (3 with A^^ = 4 and Nt = 8. All values are 
given in units of appropriate powers of the spatial lattice distance a^. 

In Ref. |1] it was shown that without optimization the Langevin dynamics 
described by Eq. fails to reproduce correct results in the limit of large 
Langevin-time. It particular, it was seen that the Langevin flow approaches 
the correct results at intermediate Langevin-times before it finally starts 
deviating. This is illustrated in Fig. [7] for the gauge invariant spatial 



20 



10 - 



' '-'plaquette''^ 
"exact"=0.91 



(Im Tr (i U|inka)/2)2 




0.1 - 



0.01 







0.1 



0.2 



0.3 



0.4 



0.5 



Langevin time 



Figure 7: The (red) solid line shows the plaquette variable Tr {7piaquette/2 defined 
in Eq. (j52p as a function of Langevin-time using g = 0.5. We employ a complex 
contour with non-zero real-time extent at temperature T = 1 as explained in the 
text. For comparison, the (black) dashed line gives the corresponding "exact" 
result for a Euclidean field theory at the same temperature. Since Tr C/piaquette/2 
is time-independent, the Langevin-time averages of both results have to agree. The 
deviation at late Langevin-times is signalled by large values for the squared link 
variable (ImTr (iC/iink i?)/2)^ defined in Eq. (f53l) . which is displayed as a (blue) 
dotted curve. 

plaquette averaged on the lattice: 



where the plaquette variable Ux,^u is defined in Eq. ([T]). The (red) solid 
line shows the result as a function of the Langevin-time for a complex 
(isoceles) triangle contour with real-time extent At^ = 1 and imginary- 
time extent Atj = 1, corresponding to a thermal theory at temperature 
For comparison, the (black) dashed line gives the "exact" value 
for this observable as obtained from stochastic quantization in Euclidean 
space time, i.e. for a time-contour with no real-time extent. In the latter 
case one can prove the convergence of the stochastic method [2j , and for the 

^See Ref. [4] for a discussion of complex time contours in this context. 




X, i<j 



(52) 



21 



employed parameters we find at late Langevin-times Tr [/piaquottc/2 = 0.91 
giving two significant digits. Of course, this comparison is only possible 
because we consider the special case of a time-independent observable, whose 
value has to be the same in Euclidean as well as Minkowskian space-time. 
From Fig. [7] we see that the result is indeed independent of the employed 
time-path for not too late Langevin-times. However, finally deviations occur 
demonstrating the breakdown of the complex Langevin method. In this 
case the plaquette variable (!52|) develops larger fluctuations of the imaginary 
part, whose Langevin-time average is zero, and in Fig. [7| we only show its 
real part. The onset time for deviations can be delayed by further decreasing 
the real-time extent of the lattice, which is analyzed in detail in Ref. [1]. In 
principle this may be used to extract physical results for sufficiently short real 
times, however, this procedure would provide severe restrictions for actual 
applications of the method. 

A characteristic measure that may be used to monitor this breakdown is 
given by the quantity 

^ (ImTr (zf/u„k o)f = ^ (ImTr (zf/,,, a)f = ^ (im^,,) , 

^ a:,i * x,j 

(53) 

where we used for the last equation the representation corresponding to 
Eq. (HHjl . This quantity is not gauge invariant and would vanish identically 
for U G SU{2). As explained in Sec. [21 for the complex Langevin dynamics 
U G SL{2, C) and, therefore, the quantity (153|1 can be non-zero and provides 
a characteristic quantity to measure deviations from SU{2). From the (blue) 
dotted line of Fig. [7]one observes that the breakdown of the complex Langevin 
dynamics occurs once this quantity becomes significantly larger than one. 

In order to make contact with the discussions of Sees. 13.11 and 13.21 we 
display in Fig. [8] the real and imaginary part of the plaquette variable 
fl52|) from snapshots with constant Langevin-time stepping. The employed 
parameters are the same as for Fig. [71 The (red) crosses give the distribution 
for sufficiently large Langevin-times, i.e. for times when the plaquette variable 
deviates from the correct results. From Fig. [71 one observes that for the 
employed parameters this is the case for -(9 > 0.5. For these times Fig. [HI 
exhibits a relatively widespread distribution of values in the complex plane. 
This is similar to what is observed in Sees. 13.11 and 13.21 for the U{1) and 
SU{2) one-plaquette models for those cases where the complex Langevin 
method fails. In contrast, for earlier Langevin-times we find a compact 
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Figure 8: Distributions of the values for the averaged spatial plaquettes 
TV f/piaquette/2 in the complex plane from snapshots of the Langevin process. The 
(red) crosses give the distribution for large Langevin-times obtained from the 
stochastic process without optimization. The (blue) stars give the distribution 
for large Langevin-times from optimized updating using gauge fixing, which yields 
correct observable averages. The employed parameters correspond to those used 
in Figs. [7] and [91 



distribution similar to what is displayed as (blue) stars in Fig. [HI Accordingly, 
for these earlier Langevin-times the values for the plaquette variable f[52l) 
agree well with accurate results. The crucial role of a compact distribution 
for the convergence to accurate results was also observed for the one-plaquette 
models in Sees. 13.11 and 13.21 In the following, we will show how to stabilize 
a compact distribution for all Langevin-times using optimized updating by 
gauge fixing. The (blue) stars in Fig. [8] actually correspond to results 
obtained from gauge-fixed Langevin dynamics, which will be explained below. 

4.2 Optimized updating by gauge fixing 

In the previous sections we have observed that the breakdown of the 
complex Langevin dynamics occurs in the presence of large fluctuations of 
the complexified dynamical variables. In the spirit of Sec. 13.21 we consider 
here a reduction of fluctuations by gauge fixing. We employ maximal axial 
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gauge, i.e. on a periodic latticqj one can fix by gauge transformations 



W-\x)U,,^W{x + fi) 



(54) 



tfie link variables to one for the following links: 



/i 



/i 



0, 0<x" < Nt 

1, 0<x^ < Ns 

2, 0<x^ < Ns 

3, < a;^ < iV, 



1, < x\x'^,x^ < N,, 
1, < x^x^ < x° = 0, 
1, 0<x'^ < Ns, x° = = 0, 
1, = x^ = x^ = 0, 



(55) 



where the lattice points are labelled by integers < x^,x^,x^ < A^^ and 
< x° < Nf. The gauge-fixed links are not updated [10] and the rest is 
updated according to Eq. (jll). 

Another possibility, which would be the analogue of what was employed 
for the one-plaquette model in Sec. 13.21 is to update all links but after each 
Langevin-time step to calculate and apply the field of gauge transformations 
W{x) in order to fix the link variables according to Eq. (!55l) . One can build 
up the field of gauge transformations by fixing first W{0, 0, 0, 0) = 1, than 
get the values 1^(0,0,0,2;) by solving {x)Ux,fj,W {x + ft) = 1 according 
to the last line of Eq. fl55l) . than the values W{0,0,y, z) according to the 
third line, and so on. The latter method turns out to be not as efficient in 
suppressing fluctuations as the maximal axial gauge. 

In Fig. [9] we show the results corresponding to Fig. Ul however, now 
obtained from the stochastic dynamics with optimized updating by gauge 
fixing. One observes that the result from the complex Langevin equation 
stays close to the "exact" result and we found no sign of increasing deviations 
at sufficiently large Langevin-times. Accordingly, also the quantity (!53l) stays 
comparably small as displayed in Fig. [91 This reflects the fact that the 
distribution of the dynamical variable in the complex plane remains relatively 
compact, which is exemplified in Fig. [S] for the plaquette variable 0521) . 

This is an important advance as compared to the results of Ref. [4J , where 
no stable physical fixed point solution of the Langevin equation could be 
observed. Gauge fixing turns out to be an efficient way to reduce fluctuations 
of the complex Langevin dynamics. However, increasing the coupling leads 
to increased fluctuations and we find that gauge fixing alone is not enough 
for sufficiently large g. Therefore, the combination of gauge fixing and not 

^The time contour is also periodic because we are studying thermal equilibrium. 
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Figure 9: Same as Fig. [7] but with gauge fixing. In contrast to the former, the 
results from the complex Langevin equation agree well to the "exact" values even 
at large Langevin-times. Note that the displayed Langevin-time exceeds the one 
of Fig. [7] by about a factor of 100. 

too large values for the coupling lead to the quantitative results we observe. 
For the results of Fig. [Hlwe use g = 0.5. For the employed parameters with 
temperature T = 1, we find that increasing g to values larger than about one 
leads to deviations from correct results at large Langevin-times similar to the 
situation displayed in Fig. [71 Decreasing the temperature or using shorter 
real-time extent improves the situation, as has been found also without gauge 
fixing in Ref. [Ij. We have checked on Ng = 4^, 32^ lattices that the value 
of g one has to use for physical results does not depend on the spatial size 
of the lattice. In principle, this means that even though the small g (and 
thus larger /3) corresponds to smaller lattice spacings, with a bigger lattice 
size one could overcome this effect at the expense of computational time. 
However, in practice this is difficult to achieve for g < 1 because of limited 
resources. 

The discussion is reminiscent of the one in Sec. 13. where it is shown that 
the (a = P) reweighted U{1) one-plaquette model is governed by classical 
dynamics in the limit /3 — ooQ The reduction of fluctuations by increasing 

^Note that the SU{2) one-plaquette model has no well defined limit /3 ^ oo as discussed 
in Sec. [321 
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(3 in that model is exemplified in Figs. and |H Qualitatively similar, for 
the gauge theory one obtains for larger (3 ~ 1/ more compact distributions 
than shown in Fig. [HI 

5 Conclusions 

For the SU (2) gauge theory with action ([2]) real-time stochastic quantization 
requires the dynamical variables to become elements of S'L(2,C). The 
expectation values of the underlying gauge theory are recovered after taking 
noise or Langevin-time averages, respectively. This change from the compact 
SU{2) gauge group to the non-compact SL{2, C) has important consequences 
for the Langevin dynamics. For the former detU = + 6^ = 1 yields 
finite values for the real a and b, using the representation of the link 
variables corresponding to Eq. (jUj). In contrast, the SL{2, C) group admits 
unbounded values for the now complex a and b. Accordingly, we observe 
that without optimized updating the complex Langevin equation yields large 
fluctuations for the dynamical variables at sufficiently late Langevin-times, 
which finally lead to a breakdown of the method. Apart from the SU{2) 
gauge theory in 3 + 1 dimensions, we observe analogous findings for the U{1) 
and SU{2) one-plaquette models. The simplicity of the one-plaquette models 
allow an analytical analysis as presented in Sees. 13.11 and 13. 2[ 

In this paper we showed that large fiuctuations of the complexified 
dynamical variables can be efficiently reduced by employing optimized 
updating procedures for the Langevin process. Here we investigated 
optimized updating using stochastic reweighting or gauge fixing, respectively. 
These procedures do not affect the underlying theory but can stabilize the 
physical fixed point of the Langevin equation. The success of stochastic 
reweighting was found to be linked to the appearance of attractive fixed 
points of the Langevin fiow, while the gauge fixing simply constrains the 
dynamical variables. For the gauge theory we demonstrated that gauge fixing 
leads to an efficient reduction of fiuctuations for not too small (3 ^ 1/ g'^: We 
employed maximal axial gauge and calculated plaquette averages on a lattice 
with non-zero real-time extent. Where applicable, the results were shown to 
accurately reproduce alternative calculations in Euclidean space-time. This 
is an important advance as compared to the results of Ref. [1], where no stable 
physical fixed point solution of the Langevin equation could be observed for 
the non-Abelian gauge theory even for small couplings. 
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For the reweighted U{1) and the "gauge- fixed" SU{2) one-plaquette mod- 
els, we obtained accurate results also for small j3. The fact that real-time 
stochastic quantization is simpler to apply to non-gauge theories is in accor- 
dance with earlier findings for scalar field theories in Ref. [4j. Since large 
P describe the continuum limit of the lattice gauge theory, with the present 
results at hand one has in principle a procedure to do simulations in non- 
Abelian gauge theory. However, without further improvements, it is difficult 
to achieve (? ^ 1 because of limited computational resources. Also further 
tests using different time contours along the lines of Ref. [3j are necessary 
before pressing questions of calculations of transport properties and nonequi- 
librium gauge theory dynamics are to be addressed. 
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6 Appendix 

In this appendix we show that the optimized updating procedure employed 
for the SU{2) one-plaquette model in Sec. 13. 2. 31 can be mapped to a Langevin 
process with one degree of freedom only, described by Eq. (ISTj) . 

Using the parameterization of Eq. PHI) we write U = a + iba. The "gauge 
fixing" condition employed in Sec. 13.2.31 then reads b = (0,0, \/l — a^). We 
denote in Eq. (H9|) the Langevin-time-stepping matrix as G = a -|- ida. 
Because U is rotated to the fixed gauge after the last Langevin-time step 
we have 6i = 62 = 0, Tr((Tit/) = Tr((j2t/) = and, accordingly, 

di = r/iv/e, d2 = r]2y/e, = -e^Ti{asU) + r/sv^ • (56) 

Multiplying and U, one obtains U' = a' + ib'a using 61 = 62 = and 

{ya){wa) = vw + ia{v x w) : 

a' + ib'a = aa — bd + iai{adi — 6/263) + ia2{ad2 + dib^) +ia3{db3 + ad^) . (57) 
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Again rotating to the fixed gauge one obtains U" = a" + ib"a with 

a" = a', b'l = b'^ = 0, h'i = ^bf + b'i + b'i . (58) 
Expanding 63 to order e according to Eq. fl56|) leads to 

One can write c?^ + = 4e + ep, where p is a noise term with zero mean. 
Since it is muhiphed with e, in the continuum hmit one can neglect this 
termj^ Using the notation 63 = sm{Lp" /2) and similarly 63 = sm{Lp/2), and 

Aip = (is + 2e^, and noticing that to order e one has cos Aip = (1 — ^) yields 

sin ^ = cos(Ay9) sin ^ + cos ^ sin( A(y9) . (60) 
Since this describes angle addition, we get: 

^ - 77 = «3 + 2e— = sm - + 2e cot - + 773 Ve , (61) 

which is the discretised Langevin equation flHTl) as obtained from the 
"reduced" action 5* = 21nsiny9/2 + if] cos(y9/2, with Langevin-time step 4e. 
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